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In the design of beam transport lines one often meets the problem of constructing a quadrupole 
lens system that will produce desired transfer matrices in both the horizontal and vertical planes. 
Nowadays this problem is typically approached with the help of computer routines, but searching 
for the numerical solution one has to remember that it is not proven yet that an arbitrary four 
by four uncoupled beam transfer matrix can be represented by using a finite number of drifts and 
quadrupoles (representation problem) and the answer to this questions is not known not only for 
more or less reahstic quadrupole field models but also for the both most commonly used approxima- 
tions of quadrupole focusing, namely thick and thin quadrupole lenses. In this paper we make a step 
forward in resolving the representation problem and, by giving an explicit solution, we prove that 
an arbitrary four by four uncoupled beam transfer matrix actually can be obtained as a product of 
a finite number of thin-lenses and drifts. 



I. INTRODUCTION 

In the design of beam transfer lines one often encoun- 
ters the problem of finding a combination of quadrupole 
lenses and field free spaces (drifts) that will produce 
particular transfer matrices in both the horizontal and 
vertical planes. Nowadays this problem is typically ap- 
proached with the help of computer routines which min- 
imize the deviations from the desired matrices as func- 
tion of the quadrupole strengths, lengths and distances 
between them. Although very sophisticated software be- 
came available for these purposes during the last decades, 
there is an important theoretical question which has 
not been answered yet and whose answer could affect 
the strategy and efficiency of numerical computations. 
Searching for a numerical solution one has to remember 
that it is not proven yet that an arbitrary four by four 
uncoupled beam transfer matrix can be represented by 
using a finite number of drifts and quadrupoles (repre- 
sentation problem) and the answer to this questions is 
not known not only for more or less realistic quadrupole 
field models but also for the both most commonly used 
approximations of quadrupole focusing, namely thick and 
thin quadrupole lenses. 

In this paper we make a step forward in resolving the 
representation problem and prove that an arbitrary four 
by four uncoupled beam transfer matrix actually can be 
obtained as a product of a finite number of thin-lenses 
and drifts. Even though our proof uses more thin-lenses 
than probably needed, we believe that the solution pro- 
vided is not only of theoretical interest, but could also 
find some practical applications because it uses explicit 
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analytical formulas connecting thin-lens parameters with 
the elements of the input beam transfer matrix. 

Though the thin-lens kick is the simplest model of 
the quadrupole focusing, its role in accelerator physics 
can hardly be overestimated. The thin-lens quadrupole 
approximation reveals the analogy between light optics 
and charged particle optics and, if one takes into account 
difficulties of analytical manipulations with the next by 
complexity thick-lens quadrupole model fl], , is an in- 
dispensable tool for understanding principles and limi- 
tations of the already available optics modules and for 
development of the new optics solutions (see, as good 
examples, papers 0-01 )■ 

The paper by itself is organized as follows. In section 
II we introduce all needed notations and give the lower 
bound on the number of drifts and lenses which are re- 
quired for a solution of the representation problem by 
providing an example of a matrix which can not be ob- 
tained using five thin-lenses and five variable drift spaces. 
This example and other of our attempts (though omit- 
ted in this paper) to find thin-lens solutions for particular 
beam transfer matrices lead us to the conjecture that in 
order to represent an arbitrary four by four uncoupled 
beam transfer matrix one needs at least six thin-lenses 
if the distances between them can be varied or at least 
seven thin-lenses if this variation is not allowed. 

In section III we prove that an arbitrary four by four 
uncoupled beam transfer matrix can be obtained as a 
product of a finite number of thin-lenses and drifts by 
giving an explicit solution of the thin-lens representation 
problem which uses equally spaced thin-lenses. The core 
idea of our approach is the reduction of the initial 2D 
problem to two independent ID problems, and the sym- 
metry of the equally spaced thin-lens system allows to 
make such a reduction while minimizing the number of 
nonessential technical details. The solution obtained uti- 
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lizes thirteen lenses if the spacing between them is fixed 
beforehand and twelve lenses if this distance can be used 
as an additional parameter. Thus it uses six more lenses 
than the minimal number stated in our conjecture, but 
the setting of these six lenses depends only on the dis- 
tance between lenses and therefore does not depend (at 
least directly) on the particular input beam transfer ma- 
trix. 

And, finally, in section IV we generalize our results 
of the previous section to the case of arbitrarily spaced 
thin-lenses. As an important by-product of this section 
we mention the representation of the matrix of the drift- 
lens system as a product of elementary P matrices (see 
formulas (pS|) - ([7T|) below). This representation could be 
a useful tool for the analytical study of the properties 
of thin-lens multiplets and it also shows the difference 
between the role of the variable drift spaces and the role 
of the variable lens strengths when they are used as fitting 
parameters. Besides that thin-lens blocks with decoupled 
transverse actions introduced in this paper are another 
point of general interest, although the idea of decoupling 
by itself is not new in the field of accelerator physics (see, 
for example, [1, Q ) . 

II. STATEMENT OF THE PROBLEM AND 
PRELIMINARY CONSIDERATIONS 

Let M be an arbitrary four by four uncoupled beam 
transfer matrix and let the two by two symplectic ma- 
trices Mx and My be its horizontal and vertical focusing 
blocks respectively. Let us denote by Q{g) the transfer 
matrix of the one dimensional thin lens of strength g and 
by D(l) the transfer matrix of the one dimensional drift 
space of length I 

Qi9) = {l',), m-{l[)- (1) 

The problem of representation of the matrix M by a thin- 
lens system can then be written as 

D{1„) Q{±gn) ■■■■■ D{h) Q(±5i) = M^^y, (2) 

where (here and later on) one has to take the upper sign 
in the combinations ± and =F together with the index x 
and the lower sign together with the index y. 

Note that the drift-lens system presented on the left 
hand side of the equations ^ consists of equal numbers 
of drifts and lenses and the first element which the beam 
sees during its passage is a thin-lens. Alternatively one 
can consider equations 

Q{±gn) D{1„) • . . . • g(±(7i) D{h) = M,^y, (3) 

where the first element is a drift space, or one can use 
the drift-lens system with a non equal number of drifts 
and lenses which starts and ends with a drift (or a lens) , 
but for the moment this is not important. 



There are many unanswered questions related to the 
equations the most interesting for us in this paper 
is the following: given a matrix M, does there exist a 
number n such that these equations have a solution? If 
the answer to this question is positive, could the number 
n be chosen independently from the input matrix M and, 
if it is also possible, what is the minimal n required? 

From a mathematical point of view the equations ^ 
are a system of eight polynomial equations in 2n un- 
knowns and for any polynomial system considered over 
an algebraically closed field of complex numbers there is 
an algorithmic way to answer the question if this system 
has infinitely many solutions or has a finite number of 
solutions, or has no solutions at all. This can be done 
by transforming the original system to a special form 
called a Grobner basis and, very loosely speaking, is an 
analogue of the Gaussian elimination process in linear 
algebra [l^. The Grobner basis can be computed in 
finitely many steps and, moreover, nowadays its calcula- 
tion can be done with the help of symbolic manipulation 
programs like Mathematica and Maple. 

Unfortunately, we are interested in the real solutions of 
the equations ^ constrained additionally by the require- 
ments for the drift lengths to be nonnegative and there- 
fore we can not use all benefits provided by the Grobner 
basis theory. Nevertheless, the use of the Grobner basis 
approach, although it did not help us to solve the prob- 
lem in general, it was very useful in providing examples 
of particular matrices which can not be obtained using 
a certain number of thin-lenses and drift spaces. For ex- 
ample, using the Grobner basis technique, it is possible 
to prove that the matrix M with 

Mx = M, = J 5 ^ (4) 

can not be represented by five thin lenses and five vari- 
able drift spaces starting either from a lens like in the 
equations ([2]) or from a drift like in the equations ([3]). 

This and many other of our attempts to study the rep- 
resentation problem for particular beam transfer matrices 
lead us to the conjecture that in order to be able to rep- 
resent an arbitrary four by four uncoupled beam transfer 
matrix one needs at least six thin-lenses if the distances 
between them can be varied or at least seven thin-lenses 
with nonzero drift spaces between them if this variation 
is not allowed. 

To finish this section, let us note that in the above dis- 
cussions we made no use of the fact that we are interested 
not in the general system of polynomial equations, but 
only in the polynomial system produced by a product of 
matrices with simple inversion properties 

Q-\9) = ^"'(0 = D{~1). (5) 

Choosing some k — 1, . . . , n — 1 and using ^ one can 
rewrite the system ^ in the equivalent form 

D{lk) Qi±gk) ■ . ■ . ■ D{h) Q{±gi) = 
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Q(Tgfc+i) D{-lk+i) • . . . • QiTgn) D{-ln) M,^y. (6) 

This trick can be used for the ehmination of a part of the 
unknowns from the original system by solving equations 
([6]) with respect to the variables gi, . . . , gkJi, ■ ■ ■ ,lk or 
one may even think to construct an iterative solution 
method which could be considered as matrix version of 
the method of successive elimination of unknowns [ll| . 
Unfortunately, however this approach did not give us any 
additional noticeable simplifications in the solution of the 
representation problem. 

III. SOLUTION OF 2D PROBLEM USING 
EQUALLY SPACED THIN LENSES 

In this section we will give an explicit solution of 
the thin-lens representation problem which uses equally 
spaced thin-lenses. Instead of equations ([2]) or ([3]) we will 
consider the system 

B(m„, ±g„, pn) ■ ■■■ ■ B{mi, ±gi, pi) = M^^y, (7) 

where as an elementary building block we take a thin-lens 
sandwiched between two drift spaces 

Bim, ±5, p) = Dip) Qi±g) D{m). (8) 

If the block length I — m + p > 0, then one can repre- 
sent the block transfer matrix in the form 

B(m, ±g, p) = S-\m, p) P{2 ± Ig) S{m, p), (9) 

where 

(-!")■ 

The definition and the properties of the matrix P (and 
other elementary matrices used in this paper) can be 
found in Appendix A. 

Let us assume that in the system ([7]) all and all pk 
are equal to each other, i.e. that 

mi = . . . = TO„ = m, pi = ...=pn^p, (11) 

and let I ~ m + p > 0. The principle simplification that 
occurs in this case is that after the substitution of the 
representation ([9]) into equations ([7]) the matrices S{m,p) 
and S~^{ra,p) cancel each other and we obtain 

P{2 ± /.9„) • . . . • P(2 ± /5i) = M,,y, (12) 

where 

Mx^y = S{m, p) M^^y S-\m, p). (13) 

One sees that while the original system ([7]) is formed 
by the product of 2n -I- 1 interleaved thin-lens and drift 
matrices (with neighboring drifts lumped together), the 



system (|12p includes only n + 2 matrices depending on 
unknowns and n of them are P matrices. 

Nevertheless, the system is still too complicated 
to find easily its solutions (or even to prove their exis- 
tence) for an arbitrary matrix AI and with the number 
of lenses n equal to six or seven as required by our conjec- 
ture. Instead we will provide an explicit solution which 
utilizes thirteen lenses if the parameters m and p are 
fixed and are independent from the input matrix M, and 
twelve lenses if m and p can be varied. The main idea 
of our solution is the reduction of the 2D problem 
to two independent or, more exactly, almost independent 
ID problems by constructing thin-lens blocks which can 
act in the horizontal and vertical planes similar to a sin- 
gle P matrix, but whose actions for both planes can be 
chosen independently. At first we will consider a solution 
of the ID problem in terms of P matrices. As the next 
step we will introduce a four-lens block with decoupled 
transverse actions and then will give an explicit solution 
of the complete 2D problem. Besides that we will discuss 
the recipe for constructing lens blocks with decoupled 
transverse actions with more than four lenses. 

Before giving the technical details let us consider one 
more example obtained with the help of the Grobner ba- 
sis technique. Let us assume that m and p are fixed and 
let the matrix M be such that the matrix M is equal to 
the symplectic unit matrix 

Mx = My = (^ liy (14) 

Then this matrix M can not be represented by less than 
seven thin-lenses and with seven lenses there are many 
solutions which geometrically can be viewed as six dis- 
tinct parallel straight lines in the seven dimensional space 
of lens strengths. 

A. ID problem in terms of P matrices 

According to our plan we will prove in this subsection 
that every real symplectic 2x2 matrix M = {rriij) can 
be represented as a product of at most four P matrices. 
First, we will consider the case of three P matrices and 
will find that three P matrices are insufficient for the 
representation of an arbitrary 2x2 symplectic matrix. 
Next we will switch to the case of four P matrices and 
will show that with four P matrices a solution can always 
be found, but it is always nonunique. 

Let us start with the case of three P matrices, i.e. from 
the equation 

P{z3) Piz2) P{zi) = M. (15) 

This matrix equation is, in fact, the system of the four 
equations for the four matrix elements and, as it is well 
known, due to symplecticity of the matrices on both sides 
of ((T5|) these four equations should be equivalent to some 
system consisting of three equations only. It is clear that 
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such a system can not be obtained simply by throwing 
away one randomly chosen equation from the system (jlSp 
and the first benefit of the special form of the matrix P 
is that the system of the four third order polynomial 
equations (|15p is equivalent to the system 



zi = -mil - m2i ■ Z3 

Z2 = -17122 

m22 ■ Z3^ -I- TO12 



(16) 



which is linear in the unknowns zi, Z2 and 23. Moreover, 
this system already has a triangular form and its solvabil- 
ity depends only on the solvability of the third equation 
with respect to the variable 23. 

Elementary analysis shows that there are three possi- 
bilities for the solutions of the system (fT6|) . If 11122 7^ 0, 
then there exists a unique solution 



Zl 



17121 - 1 



Z2 



-m22, Z3 



TO12 + 1 

^122 



(17) 



If m22 = and 77121 = 1 (i.e if M = — F(— mn)), then 
there exists a one-parameter family of solutions 



-Zl + 23 = -TOll, 22=0. 



(18) 



And, finally, if TO22 = and r7i2i 7^ 1, then there is no 
solution at all. 

Very loosely speaking, the condition m22 — defines 
the two dimensional surface of singularities in the three 
dimensional space of 2 x 2 real symplectic matrices. This 
surface, in the next turn, contains the one dimensional 
curve selected by the additional relation m2i = 1. If the 
matrix M (represented as a point in our three dimen- 
sional space) lies outside of the surface of singularities, 
then a solution for such a matrix exists and is unique. 
If the point representing the matrix M belongs to the 
surface of singularities, then we either have many solu- 
tions or none depending on whether this point lies on the 
above defined one dimensional curve or not. 

Let us now turn our attention to the equation 



F(24)P(23)F(22)P(2l) = M, 



(19) 



which includes four P matrices. The equivalent to this 
equation system is given below 



21 = m2i - (mil + m2i • 24) • 23 

22 = -TO12 - TO22 ■ Zi 

(mi2 + 77122 • Z4) ■ 23 = 77122 - 1 



(20) 



and the easiest way to obtain it is to substitute into the 
system (fT6|) the elements of the matrix ^^^(24) • M in- 
stead of the TTiij . 

The system (|20p is not linear anymore, but still has a 
triangular form and its solvability depends again only on 
the solvability of the third equation with respect to the 
variables 23 and 24. Because the matrix M is nondegen- 
erated its elements 77112 a-nd TO22 can not be equal to zero 
simultaneously and therefore the expression mi2-l-r?i22-24 



considered as a function of 24 can not be equal to zero 
in more than one point. It means that the last equation 
in ([20I) always has solutions and a good way to under- 
stand their complete structure is to consider this equa- 
tion as the equation of a curve on the plane (23,24). If 
77122 • {i7i22 ^ 1) 7^ this curvc is a hyperbola with two sep- 
arate branches, if 77122 = 1 it is a degenerate hyperbola 
consisting of two intersecting lines 23 = and 24 = — 77112, 
and, finally, if 7*722 = we have a single straight line 
23 — —m^2- So we see that with the help of the four P 
matrices a solution of our problem can always be found 
and is always nonunique. 



B. Four- lens block with decoupled transverse 
actions 

Let us denote by W^ ^ the following combination of 
four P matrices 

W^'y = P{2 ± lgi)P{2 ± lg3)P{2 ± lg2)P{2 ± ^51), (21) 

which in the original variables ([7]) includes four thin- 
lenses (four-lens block). 

If one chooses (5 = ±1 and if one takes 
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93 



1 



then the block matrix can be written as 
where A is a diagonal scaling matrix with 

u^'^y ^ 2t sVi, -uy = 

and 

= 7 + • Igi + • Igi 
vjy = 7 — ■ Igi — uy ■ lg4 



Since for any given value of and to*' the equations (|25|) 
and ([26]) can be solved with respect to the variables gi 
and 34 



(22) 
(23) 

(24) 

(25) 
(26) 



9i 



54 



S uy ■w'' + u"" -wy - 28 

5 u'' ■w'' + uy -wy - 28 



I 



8V3 



(27) 



(28) 



the formula (1^5)) gives the result which we were looking 
for. Both matrices and Wy are similar to a single 
P matrix (with an inessential minus sign) and both pa- 
rameters and vjy can be chosen independently by an 
appropriate setting of the first and the last lenses in the 
block. 
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C. Reduction of 2D problem to two independent or 
almost independent ID problems 

Since with four P matrices we always can solve the 
ID problem, let us first consider a combination of four 
blocks of the type Using (|JT7)) one can show that 

the total matrix of this sixteen lens system can be written 
as follows 

Wt" W^'^ W^'^ W^'y - A (a"'^) • 



at the solution announced in the introduction, namely 
thirteen lenses if the spacing between them is fixed and 
twelve lenses if this distance can be used as an additional 
parameter. Below we will consider in detail the case of 
twelve lenses (three blocks) with variable spacing and the 
check that the use of an additional lens for the fixed spac- 
ing also works we leave as an exercise for the interested 
reader. 

In analogy with (|29|) the combination of three blocks 
can be written as 



P iw^'^) P {w^s'^) P {w^'^) P (wl'^) A (a"'^) , (29) 



w^'y w^'y wi^y 



where 



A-i (a"'2/) p p p ^ (34) 



and 



x^y x,y 
„x,y 



'■1 "3 



1*2 U4 



(30) 



■,y x,y ' ^2 ' I (31) 



where 



and 



/ x^y x,y 
^x.y _ / ^1 ^3 



(35) 



n^^y — 



^3 ' 



x.y 



(32) 



w 



x,y "'2 x,y 



(36) 



Plugging this representation into equations we ob- 
tain 



Wo 



x,y 
x.y 



,x,y 



(37) 



p [wi^y) p {wi^y) p [wi^y) p [wi^y) 



(33) 



Let us choose arbitrary nonnegative m and p with / ~ 
m + p > Q and select for each four-lens block its own 
5 — ±.\. This, in accordance with the formulas ([22]) . gives 
us the setting of the eight lenses in our system and this 
completely determines the matrix on the right hand side 
of the equations ([33|). As the last step we take 1x1% and 
as some solutions of two independent ID problems of 



the type (|19p and define the strengths of the remaining 
eight lenses using the formulas ([31])- ([32]) and ((27)) - ((28)) . 

One sees that using four blocks with decoupled trans- 
verse actions the complete 2D problem can always be 
reduced to two easily solvable independent ID problems. 
But do we really need four blocks for making such a re- 
duction? The answer is no and the reason for this is as 
follows. We know that for most of the 2x2 symplectic 
matrices the ID problem can be solved with three P ma- 
trices, which means that for most of the 4x4 uncoupled 
beam transfer matrices the 2D problem can also be solved 
with three blocks. The problem is what to do with the 
rest? Happily it turns out that by appropriate choice of 
the parameters m and p one can always move the input 
matrix M away from the region of unsolvability and, if 
the variation of m and p is not allowed, this can be done 
by using only one additional thin-lens. Thus we arrive 



x,y 

h 

x.y 

U2 



(38) 



Plugging again this representation into system ((12)) we 
obtain the equations 



P iw^^'«) P (w;'^) P iw'l^'') 



(39) 



We know that the sufficient condition for these equa- 
tions to be solvable with respect to the unknowns w^'*' 
is that the horizontal and vertical parts of the matrix on 
the right hand side both have nonvanishing r22 elements. 
The direct calculation gives us 



x.y 

mmii 



x.y 
■pm^2 



x.y 

mp 17121 



m + p 



(40) 



where m^^'^ are the elements of the input matrix M. 

Looking for a solution one can proceed further in the 
same manner as in the four block case with only one 
difference. At the first step one has to take not arbitrary 
nonnegative m and p, but such m and p that both rfj 
and are nonzero, which due to symplecticity of the 
matrices and My is always possible. 
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D. Recipe of construction of lens blocks with 
decoupled transverse actions 

In this subsection we give the recipe for the construc- 
tion of lens blocks with decoupled transverse actions. As 
we will see, this recipe works not only for the four-lens 
combination considered above, but is also applicable to 
blocks with a larger number of lenses. 

Let us consider g— lens block with g > 4 



P(2±/5,).....P(2±/5i), 



and let us assume that the product of the {q 
matrices in our block takes the form 







(41) 
2) inner 

(42) 



Then, as one can show by direct multiplication, the total 
block matrix becomes 



= -sign(-u^'^)- 



A- 



where 

„x,y 2 ± Ig 



w 



' [VW^) P{w-^y)^ (/m) , (43) 
- + lu^'^l (2 ± Igq) + sign(M"^^) e"'^. (44) 



If for arbitrary given values of and the equations 
can be solved with respect to the variables gi and gq, 
then it will be exactly what we need, and the necessary 
and sufficient condition for such solvability is 



^ \uy 



(45) 



So, in order to construct the g-lens block with the decou- 
pled transverse actions, one has to solve two equations 
making the rn elements of the x and y parts of the prod- 
uct of the {q — 2) inner matrices equal to zero and one 
has to satisfy one additional inequality constraint (j45p . 

The solution for the four-lens block was already given 
above and is unique up to a sign change (5 — ±1). Let 
us now consider more complicated (but still analytically 
solvable) case of five lenses. In this situation all possible 
solutions which bring the product of the three inner P 
matrices 



P(2 ± Ig^) P{2 ± Igs) P{2 ± Ig^) 



(46) 



to the form (|42p can be expressed as a function of pa- 
rameters I and 53 as follows: 
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1 Igz + 6^{{lgz? - 2) ■ {{2lg^f - 9) 

I ■ {ig^f - 3 



(47) 



1 ^.93 - V((^.93)^ - 2) ■ {{2lg^Y - 9 



(5 = ±1, and I > and 173 are such that 



I93 e (- 



-00, -%/3 ) U ( -\/3, -1.5 



U 



V2, %/2 U 1.5, VS) U (V3, +00 



(49) 



To complete the block construction we have to select from 
all these solutions a subset on which the functions 



1 - {lg2 T 2) • (/53 + Igi 



(50) 



satisfy the inequality (|45|) . As one can check, this can 
be achieved simply by removing from the set (j49p the 
endpoints of the given set intervals, i.e. by removing the 
points ±1.5 and ±\/2. So we see that there are many 
solutions which allow us to construct from five lenses the 
block with decoupled transverse actions and for selecting 
one of them some additional optimization criteria could 
be involved. 

Note that in the blocks constructed according to our 
recipe the setting of the internal lenses does not depend 
on the setting of the first and the last lenses and depends 
only on the geometrical block parameters (distances be- 
tween the lenses), which will be seen more clearly in the 
following section where we will consider the case of arbi- 
trarily spaced thin lenses. 

Note also that the horizontal and the vertical matri- 
ces between the first and the last lenses, when calculated 
using the original representation ([7]), both have ri2 ele- 
ments equal to zero (i.e. the phase advances between the 
first and the last lenses in the block are always multi- 
ples of 180°), but this alone without the inequality (05]) 
satisfied does not give us the block with the decoupled 
transverse actions. 



IV. GENERALIZATION TO THE CASE OF 
ARBITRARILY SPACED THIN LENSES 

When the distances between the lenses are not equal 
to each other, we immediately lose the advantage of the 
cancellation of 5* matrices between the P matrices after 
substitution of the representation © into equations ([7]) . 
Nevertheless, as we will show below, this case can also be 
treated with the tools developed in the previous section. 

Let us denote by dfci.fej the distance between the lenses 
with the indices ki and ^2 (fci < ^2). We start from the 
observation that for k = 2, . . . ,n the following identity 
holds 



S{mk,Pk)S '^{irik-i, Pk-i) 



L 



Ik 



-1 A 



lk-l,k 



dk-l,k 



U 1- 



Ik-l.k 



(51) 



which can be shown by direct multiplication and which 
requires that all Ik and dk-i,k are positive. 
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Let us now substitute the representation ([9]) into the 
equations (O and then plug in in the corresponding 
places the right hand side of the identity ([51]) . After 
that the property (|A19I) allows us to eliminate from the 
result all L and U matrices while shifting their argu- 
ments to the arguments of the neighboring P matrices, 
and leaving us with a product consisting of alternating 
P and A matrices. Although the A matrices can not be 
eliminated completely, they can be moved either on the 
left or on the right hand side of all P matrices with the 
help of the property (IA17P . As the last step we transfer 
all matrices from the left and right sides of the obtained 
solid block of the P matrices to the right hand side of 
our equation, hide them in the matrix Mx^y and end up 
with the equations 



p{vi'y)-...-p{vi^y) 



Mx 



(52) 



which already have the desired form. The detailed struc- 
ture of the arguments w^'^ and of the matrix Mx^y de- 
pends on the particular ways how the individual A ma- 
trices were moved (to the left or to the right sides) and 
is given below for the case when during transformations 
all A matrices were moved to the left hand side of the P 
matrix block. Nevertheless, the expressions given below 
are general in the sense that they contain an arbitrary 
positive parameter ci , and with the proper choice of this 
parameter one can account for all possible ways of move- 
ment of the individual A matrices. 

Mx^y = A(c„) S{mn,pn) Mx,yS~'^{mi,pi)A{ci), (53) 



dk-i,k 



{q — 2) inner matrices to the form (l42l) while also sat- 
isfying the inequality constraint (|45p . For the four- lens 
case 

w^^y = p{i'i^y)p{ii'y)p{ii'y)p{vi^y) (58) 

the two equations making the r\\ elements of the x and 
y parts of the product of the two inner matrices equal to 
zero are 



1, 



and have a solution 
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dia V 0^2,3 d2A ' 




(59) 



(60) 



(61) 



which again is unique up to a sign change [5 — ±1). The 
values u^^y for this solution are 



-'3,4 



^T^-J^-^ .(62) 

"2,3 V «2,3 ai,3 / 



Both of them are positive and clearly satisfy the inequal- 
ity (|45|) . With this choice for g-2 and gz the total block 
matrix takes the form 



where 



'(u;^'^) A (V^) (63) 



,x,y . 



1, 



(64) 



dk-i.k 1 



^/h-i Ik Ck-i ■ 



fc = 2,...,n, (55) 



ci is an arbitrary positive parameter and, because we 
do not have lenses with indices and n + 1, we use the 
conventions that 



d(3,i — h, do^2 — do^i + di^2, 



(56) 



and all results of the previous section concerning the re- 
duction of the 2D problem to two ID problems become 
applicable with some minor changes connected with the 
difference in the matrices Mx,y and Mx,y defined by the 
relations ([T3l) and (|53p respectively. Note that if, when 
placed in the beam line, the actual decoupling block 
starts from the lens with the index k, one has simply to 
add fc— 1 to the indices 1, 2, 3 and 4 in all above formulas. 



dn,n+l — In, C^n-l.ri+l — f^n- l.n + f^n.n+l • (57) 

Note that for Ci = 1 and with the condition ([TT]) satisfied 

the representation ([5^ turns into the representation 

as one can expect. Note also that by proper choice of ci 

one can for even n satisfy c„ — ci and for odd n satisfy 
_ -1 

Now in order to continue we need a lens block with 
the decoupled transverse actions and, as it is not diffi- 
cult to check, the recipe given in the previous section is 
applicable without any changes. For the construction of 
the g-lens block we still need to bring the product of the 



A. Removing of superfluous parameters 

The equations (O contain 2n parameters which specify 
the drift lengths (mi,pi, . . . , to„,p„) while only n -|- 1 
parameters, namely ?Tii,di.2, ■ . . ,dn-i,mPn have a clear 
physical meaning and are independent. Let us have a 
closer look at the formulas ([52|) - (f57| and count how many 
superfluous parameters are still left in them and then 
show ways how to remove them. 

The superfluous parameters pi and rn„ are clearly 
present, either directly as the arguments of S matrices 
or through the lengths of the first and the last building 



8 



blocks li and Z„. And actually that is all. The presence 
of the other superfluous parameters through the values 
I2, ■ ■ ■ , In-i is completely imaginary. To show this let us 
note that these values can enter the main formulas ([5^ - 
(|54p only through the values ci and c„ and through the 
combinations cf Zi, . . . , c^ln- So if we choose ci to be in- 
dependent from ^2, • • • , In-i, then these parameters can 
enter in none of the combinations c^fc due to the recur- 
sion relation 

clh = dl_^i^- fc = 2,...,n, (65) 

which follows from the recursion relation ((55|) . and like- 
wise they can not enter the value c„ because one can 
write that c„ = y^c^ln/ln- 

Thus there are only two superfluous parameters, pi and 
TO„, present in our formulas, either directly or through 
the values h and In- Do we need to remove them? In 
general not, because it is clear that none of the phys- 
ically meaningful answers will depend on them and, in 
this sense, their absence in the final results (like in for- 
mulas (|60| and (|6T|)) could work as some indirect indi- 
cator of the correctness of the calculations. But from 
another point of view, it seems better not to have any 
superfluous parameters from which one can expect noth- 
ing except some possible additional complications. 

The simplest way to remove the parameters pi and to„ 
from the formulas ((52)) -([54 | is to make them functions of 
the physically meaningful parameters. For example, one 
can take pi = 0.5 • di^2 and m„ = 0.5 • (i„_i^„. How- 
ever the way which we prefer is the modification of the 
formulas ((52|) - (l54l) in such a way that the superfiuous pa- 
rameters will disappear automatically. In doing so let us 
first present the final result and then make some remarks 
on how it can be obtained. 

P«'^).....P(z;r^) = M.,„ (66) 
Mx,y = Jk~\hn)U{-pn)M^^yU{-mi)K{bi), (67) 

^A- *= 2.-^."- 1.(69) 

\Uk-l,k Uk,k+1 J 

v:^y^bl(-^±g,y (70) 

bi > 0, bk = dk-i,k ■ fc = 2, . . . , n. (71) 

bk-i 

In order to obtain these formulas from ([52|) - ([57l) let us 
first introduce the parameters b^. = Ck^/Tk and then as- 
sume that ci is chosen in such a way that 61 does not 



depend on any superfluous parameter (for example, one 
simply can take Ci — 1/ \/Ti) ■ After this one sees that the 
parameters li and Z„ enter the left hand side of the equa- 
tion ([5^ only through the matrices P(?;^'^) and P('D^^'^). 
Due to the property (|A19p these matrices can be decom- 
posed into the following products 

Pdl'V) = P(vl-y) L(cl) = P(vl'y) L(b\ I h), (72) 

P(«i^'^) = U{-cl) PK'y) = U{-bl I Q P{v^^y). (73) 

As the last step, one has to substitute these decomposi- 
tions back into the equation ([5^ . transfer U and L to the 
right hand side and, after some straightforward manipu- 
lations, arrive at the final result described in the above 
formulas ((66)) - ([7T|) . 

Note that the whole story about the presence of the 
superfluous parameters is the result of our desire to have 
the expressions for the problem description (expressions 
([5 ^ - ([57)) ) which reduces to the highly symmetric expres- 
sions ([T^ and (IT^ in the limit of equal distances between 
thin lenses. If one does not require that, then, as we will 
outline below, it is possible to arrive at the representation 
dnHl-lini) without using the identity ©. 

According to (|A20|) and (jA21l) the matrix of the build- 
ing block can be written as 

P(m, ±g,p) = P{-p) P{±g) P{-m) J. (74) 

Substituting this representation in the original equations 
© and using that due to (IM)) 

P(-mfc) JP(-pfc_i) = -P(-dfc_i,fc) (75) 

we obtain 

P{±gn)P{-dn-l,n) ■■■■■ P(-di,2)P(±.9l)A(6i) = 

(-1)"- VC/(-p„) M^^y Ui-mi) A(6i), (76) 

where we have already introduced an arbitrary positive 
parameter bi. Now, assuming that all distances between 
lenses are positive and using (|A16I) . we can replace for 
each k = 2, . . . ,ji the matrix P{—dk-i.k) by the matrix 
—K[dk-i^k) with simultaneous adding to the arguments 
of the two neighboring P matrices the value d'j^^^ j,. After 
these manipulations we arrive at the expression 

-PKT-l,n ± .9n) A(d„_i^„)- 
P{dn-l^n + dn-2,n-l =•= 9n-l) A(c?„-2,n-l) ' ■ ■ ■ 

.... P(d2-j + ± g2)K{d^a)P(d^\ ± .9i)A(6i) = 

JU{-pn)M,,yU{-m^)k{bi), (77) 

and the last step, which is still necessary in order to ob- 
tain formulas (I551) - (|7ip . is to move all A matrices to the 
left in the left hand side of the equation ((77|) using the 
identity (jA17l) with a subsequent transfer of the matrix 
A(&~^) from the left to the right hand side of the obtained 
equality. 
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Appendix A: Elementary matrices and their 
properties 

The elementary symplectic P matrix which is defined 
as follows 



P(a) 



a 1 

-1 



(Al) 



and which we use extensively throughout this paper was 
found empirically by the usual trial and error method 
during attempts to reduce the problem of analytical 
study of thin-lens multiplets to some "more manageable" 
form. As we will see below, this matrix possesses many 
interesting properties not only by itself, but also in com- 
bination with the other elementary matrices. Although 
not widely known in the scientific community, it was no 
surprise, as we found later, that it was successfully used 
in some special area of abstract algebra ,12]. 

In order to give an expression for the product of n 
elementary P matrices let us first define a sequence of 
polynomials k„ in the variables zi , . . . , z„ recursively by 
the following equations 



K_l =0, kq = 1, 



(A2) 



I (-21, ■ 



Kn-2{zi,. . . ,Zn-2), n > 1. 

With these notations we assert that 
P(a„)-...-P(ai) = 



(A3) 



K„(ai, . 
-K„_i(ai, 



■ ,a„) 



Kn-l(a2, • ■ • , a„) 
l) ~Kn-2(a2, • ■ • , 



,(A4) 



which is clear for n = 1 and in general case can be proven 
by induction. Because such induction can be made in two 
different ways, either by adding one more P matrix from 
the left or from the right side, it is easy to see that the 
polynomials also be defined by (jA2[) and by the 

recursion relation 

{zi , ■ ■ ■ , Zfi) = Zi ■ Kn~l{z2, ■ ■ ■ t Zn) 



Kn~2{z3, ■ ■ ■,Zn), n > 1. (A5) 

Comparison of (jA3p and (|A5p implies that 



According to (|A4p we can write down the matrix of the 
product of any number of elementary P matrices without 
making any matrix multiplications. In this connection 
let us note that the problem of deriving some recursion 
relations which allow to obtain the transfer matrix of an 
arbitrary multiplet without actual matrix multiplications 
was also addressed in T]. 

It is clear that the matrix P(0) coincides with the 2x2 
symplectic unit matrix J, i.e. that 



P(0) 



1 

-1 



J, 



(A7) 



and the following relations between the P matrices can 
be easily verified by direct multiplication: 



P{a)JP{b) = -P{a + b), 



P3(±l) = Tl, 



(A8) 



(A9) 



p-i(a) = JP{-a) J ^ a-I - P{a), (AlO) 



P{a)p-\b) = -P{a-b)J, (All) 



P{a) P~^{b) P{c) ^ P{a-b + c), (A12) 

where / is the 2x2 identity matrix. 

Let us now introduce three more elementary matrices. 
The diagonal (scaling) matrix 



A(«) 



a 




(A13) 



and the lower and upper triangular matrices with unit 
diagonal elements 



L{a) 



1 
a 1 



U{a) 



1 a 
1 



(A14) 



Note that although the matrices L and U formally coin- 
cide with the matrices of the thin lens and the drift space 
respectively, we have introduced them in order to distin- 
guish the situations where matrix of lens or drift has 
physical meaning and where the usage of low or upper 
triangular matrix is simply the reflection of the mathe- 
matical technique used. 

We have the following relations between the matrices 
P, A, L and U: 



P{a) P{a~^) P{a) = A(-a), 



(A15) 



P{a)P{b-^)P{c) = P{a-b)K{b-^)P{c~b), (A16) 



A(a)P(&)A(a) = P{aH), 



(A17) 
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P{a-^) = L{-a)A{a-^)U{a), (A18) U{a) = -P{-a) J. (A21) 

U{a) P{h) L{c) = P{b + c- a), (A19) 

Although these relations are elementary, they are basic 
L{a) = —JP{a), (A20) for all results of this paper. 
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